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Abstract 

The calculation of correlation functions for j3 = 1 random matrix ensembles, which can be carried 
out using Pfaffians, has the peculiar feature of requiring a separate calculation depending on the 
parity of the matrix size N. This same complication is present in the calculation of the correlations 
for the Ginibre Orthogonal Ensemble of real Gaussian matrices. In fact the methods used to compute 
the P = 1, N odd, correlations break down in the case of N odd real Ginibre matrices, necessitating 
a new approach to both problems. The new approach taken in this work is to deduce the j3 = 1, N 
odd correlations as limiting cases of their N even counterparts, when one of the particles is removed 
towards infinity. This method is shown to yield the correlations for N odd real Gaussian matrices. 
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1 Introduction 



1.1 Background 

The identification of the statistics of a complicated physical system with the statistics of random matrices 
was first made explicit by Eugene Wigner [29j SO] ■ Motivated particularly by the difficulty of applying 
an individual particle model to the calculation of nuclear energy levels [16 , he studied the eigenvalues 
and eigenvectors of an ensemble of real symmetric matrices whose entries were normally distributed. 

Due to concerns with defining a unique Gaussian ensemble and the difficulty of carrying out the 
mathematics, Dyson [Jj studied Circular Ensembles, where the matrix elements are points on the unit 
circle. In Dyson's formulation, he hypothesised identifying the behaviour of a sequence of n eigenvalues 
(from a total of N 3> n eigenvalues, all lying on the unit circle), with the behaviour of n energy levels 
in a physical system. As part of the analysis, he developed the seminal distinction between Orthogo- 
nal, Unitary and Symplectic ensembles by examining time-reversal and rotational symmetries [B] — the 
required symmetries imply that the respective ensemble be invariant under orthogonal (COE), unitary 
(CUE) and symplectic (CSE) transformations. 

The analogous ensembles in the Gaussian case are the GOE, GUE and GSE. The required symmetries 
mean that GOE consists of real symmetric matrices, GUE of Hermitian complex matrices, and GSE of 
self-dual real quaternion matrices. As a consequence of these facts the eigenvalues of all the Gaussian 
ensembles are constrained to lie on the real line; the joint probability density functions (jpdfs) are [23] 

P Nt ^ Xl ,... )XN ) = J— e-^=^^ 2 T\\ Xi - Xj f (l) 

with (3 — 1,2,4 corresponding respectively to the GOE, GUE and GSE; Cn,p is some normalisation 
constant. 

Compare this with the jpdf for the circular ensembles: 

Q N ,p{6 u ...,6 N ) = -L_ J[ | e ». _ f ( 2 ) 

where (3 = 1,2,4 correspond respectively to COE, CSE and CSE. Cup is again used to denote the 
normalisation. The eigenvalues of these ensembles lie on the unit circle in the complex plane. Note 
that (fTJ) and ([2]) share the common feature that the interaction between the eigenvalues is precisely the 
distance between them raised to the power of (3. 

The clear distinction between the elements of the matrices of each Gaussian ensemble led Ginibre [14] 
to a generalisation specified by removing the symmetric/Hermitian/self-dual requirement. Despite the 
irrelevance of the orthogonal, unitary and symplectic terms in relation to their formation, these ensembles 
are named for their eponymous Gaussian cousins: GinOE of real matrices; GinUE of complex matrices; 
and GinSE of real quaternion matrices. 

Ginibre's ensembles provided the motivation for the current work. GinUE is relatively straightforward; 
Ginibre himself worked out the the jpdf and eigenvalue correlations in his original paper of 1965. GinSE 
was slightly more opaque; Ginibre was able to identify the jpdf but was not able to find the correlations. 
Further work has since developed these correlations [15] . On the other hand, GinOE proved quite 
intractable until very recently. Ginibre was only able to find the jpdf in the case where all eigenvalues 
are real. The full jpdf was not developed until 1991 [17] . and the complete correlations took longer still 

m- 

The trouble largely stems from there being different sectors to the jpdf, determined by the number of 
real eigenvalues [17j HI] . Yet there is a further problem. 

For an TV x N matrix, the N even and ./V odd cases require a different analysis. In the N even case 
the eigenvalue correlations have been the subject of the recent works [3] [TT] [5S] [57] . Of these, only [57J 
gives an analysis of the correlations for N odd. The feature that N odd and N even must be treated 
separately is already present in the calculation of the correlations for Jl} and ^ in the case (3 — 1. It is 
the aim of this paper to show that the calculation of these latter correlations, and those of GinOE, can 
be accomplished as a limit of the N even cases. 
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1.2 Parity problems for (3 = 1 

For (3 — 1 in (|TJ) or @ the general n-point correlation is given by an n x n quaternion determinant, or 
equivalently a 2n x 2n Pfaffian (these are revised in Section [2]) . This fact was revealed by Dyson [7J in 
the case of © and further developed in a series of papers by Mehta [30] [3U HH] • The strategy of Dyson 
was to express the jpdf itself as an N x TV quaternion determinant, and then to successively integrate this 
to deduce the n-point correlation. Successive integrations require a certain skew-orthogonality property 
that effectively reduces the size of the quaternion upon each integration. It is the first step which requires 
that N even and N odd be treated separately. In the method of Dyson the identity between the jpdf and 
the quaternion determinant involves a matrix factorisation; in the N odd case this was achieved by the 
introduction of an auxiliary parameter 5, and the limit 6 — > is taken. 

These complications make N odd technically more difficult. In fact, the study of the correlations for 
a generalisation of (JTJ with 0=1 

±f[e-y™ I] (3) 

N 2 = 1 l<j<k<N 

for general V(x), given in [18j . was restricted to N even. The modification required for N odd is implied 
in the later paper of Frahm and Pichard [13j , and the resulting formulae were written out explicitly in 

The method of integration over alternate variables [19 , or use of integration formulae due to dc Bruijn 
for the product of a Pfaffian times a determinant [5] (each of which implies distinct treatment of even 
and odd), can be used to express the generating function for the correlations as an N x N Pfaffian. In 
the TV even case de Bruijn's method was used by Tracy and Widom [28] to obtain a formula for the 
generating function as the square root of the Fredholm determinant of a 2 x 2 matrix integral operator, 
and functional differentiation is used to extract the correlations. Using essentially this same method 
Borodin and Sinclair computed the correlations for GinOE in the case of N even. 

For the case of N odd, a quite different method, using Grassmannians and without the use of skew- 
orthogonal polynomials, was developed in [27) . By introducing an artificial Grassmannian in the odd 
case, and using a diagrammatic method of expansion, the nth order GinOE correlations are generated 
for general N. 

1.3 Guide to paper 

Section 2 introduces the terminology and concepts required for the study of /3 = 1 cases, and outlines 
the problems encountered when dealing with an odd number of eigenvalues. The difficulties that have 
postponed the identification of the jpdf and the calculation of the correlation functions for the GinOE 
are described, along with some of the attempts to overcome them. 

Section 3 details the process of developing N odd correlations for the jpdf (J3j) as limiting cases of the 
extant N even correlations pQ. See also [10]. 

In Section 4, the method as presented in Section 3 is applied to the GinOE using the existing N even 
solution. The final section discusses the conditions under which this method would be applicable and 
suggests a future application. 

2 The trouble with being odd 

Dyson [7J identified that the correlation functions for j3 = 1,4 are naturally expressed as quaternion 
determinants of matrices with quaternion elements. The problem can also be looked at from a Pfaffian 
viewpoint, since they are equivalent. This is the essence of the parity problems; Pfaffians of odd dimension 
are trivially zero, and their machinery does not seem immediately adapted to dealing with this. 
We first look at the relevant definitions. 
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2.1 Quaternion determinants and Pfaffians 

For a quaternion expressed as a 2 x 2 matrix 



a b 
c d 



the dual of M iy j is written as 

d —b 



(4) 
(5) 



-c a 

A matrix M consisting of TV x TV blocks of quaternions is said to be self-dual if 

Mj,i = M id (6) 
For ease of reference we provide here the definitions of quaternion determinants and Pfaffians. 

Definition 2.1. Quaternion determinants 

Let M be an N x iV self dual matrix of quaternions M i;J - . The quaternion determinant is defined by 

l 

QDet[M] = ]T (-l) N - l l[(M ab M bc ...M sa )W ( 7 ) 

Pes N l 

The superscript (0) denotes the operation ^Tr of the quantity in brackets. P is any permutation of 
(1, ...,N) which consist of I disjoint cycles of the form (a — > b — > c — > • • • — > s — > a). If the Mjj are scalars 
then QDet[M] = Det[M]. 

Definition 2.2. Pfaffians 

Let X = [Tij]ij—i i ... i 2N where Tji — —nj so that X is an antisymmetric matrix of even degree. Then 
the Pfaffian of X is defined by 

Pf[X] = £(-P) T P(l),P(2) T P(3),P(4) ' ' ' T P(2AT-1),P(2AT) 

P(2i)>P(2/-l) 

= 2^V! 51 £ ( P ) T P(1),P(2)TP(3),P(4) ' ' ' r P{2 N-l),P(2N) (8) 

' Pes 2N 

where e(P) is the sign of the permutation. The * above the first sum indicates that the sum is over 
distinct terms only (that is, all permutations of the pairs of indices are regarded as identical). 

Usefully, Pfaffians can be calculated using a form of Laplace expansion. For a determinant, recall 
that you may expand along any row or column. For example, expand a matrix A = [ay]i 1 _ 7 -=i ) ... n along 
the first row: 

Det[A] = OLiDetlA] 1 - 1 - ai >2 Det[4] 1,2 + • • ■(-l) n+1 oi, n Det[i4] 1>n (9) 

where Det^] 1 '- 7 means the determinant of the matrix left over after deleting the ith row and jth column. 

The analogous expansion for a Pfaffian involves deleting two rows and two columns each time. For 
example, expanding a skew-symmetric matrix B = [6y]i 1 _ 7 -=i ) ...„ (n even) along the first row: 

Pf[S] = ftLiPffB 1 ' 1 ] - 6i, 2 Pf[B 1>2 ] + • • ■(-l) n 6i,„Pf[B 1>n ] (10) 

where Pf[73 JJ ] means the Pfaffian of the matrix left after deleting the ith and jth rows and the ith and 
jth columns. Laplace expansion requires n\ calculations for a determinant, and n\\ = n- (n — 2) • (n — 4) • •• 
in the case of a Pfaffian. 
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2.2 A Pfaffian by any other name 



Define the matrix 



Z- 



2N 







(11) 



(that is, Z2N is a 2N x 2N matrix consisting of N 2 x 2 anti-symmetric blocks along the diagonal). By 
operating with Z2N we can change an N x N self-dual quaternion matrix into a 27V x 2N anti-symmetric 
matrix. The effect of Z-?n is to interchange each pair of columns and multiply every second column by 
(-!)• ^ " 

Using Z 2 n we find the following simple relationship between Pfafhans and quaternion determinants 



QDet[M] = Pf [MZ. 



2NI 



(12) 



For a proof of {T2J) see [10] . 

In light of Q12[) it is understood that we may use the terms quaternion determinant and Pfaffian 
interchangeably. 



2.3 2N or not 2N 

Why is N odd for (3 = 1 special? This can best be understood by considering the generating function 



J N,0[a\ 



N 



dx N Y[a(xi)e-^^ J] \x k - Xj f 

l<j<k<N 



1 = 1 



In the case /3 = 4 the generating function can be expressed in the Pfaffian form [18] 

Z NA [a] = Nl 2 N Ff[a jtk ] jtk=ll 



,2N 



where 



1 



e- 2y ^a(x)(Q,_ 1 (a;)g fc _ 1 (x) - Q fc _ 1 (a:)O i _i(a!))da! 



with Q n (x) an arbitrary monic polynomial of degree n and 

Qn-i{x) := ^(e- 2V ^Q n -i(x) 
In the case (3—1 and TV even, the method of integration over alternate variables gives [18] 

Z NA [a] = Nl 2 N ' 2 Pf[7i,fe]j, fc =i,...,iv 

where 

1 

7i,fc ~ 2 



e-^^a^fl^i^) / e-^^a^)^!^) sgn(y - as) 



(13) 



(14) 



(15) 



(16) 



(17) 



(18) 



with Rn(x) an arbitrary monic polynomial of degree n. The structure of (|14| is subtly, but crucially, 
different to (|17j) as the matrix in the latter is of even size only when N is even, while that of the former 
is always of even size, regardless of the parity of N. 

According to Definition 12.21 a Pfaffian is defined only for an even sized matrix. Extension of this 
definition to the odd-sized case has been given by de Bruijn [5] — his treatment amounts to bordering 
by an additional column of Is and an additional row of —Is with a zero in the bottom right corner. In 
fact the method of integration over alternate variables also gives that the N odd Pfaffian formula [22] is 
obtained from the one for N even (|17p by bordering 



2 Nll [a]=m 2< Ar + 1 )/ 2 Pf 



[yi,j]i,j=i,...,N [^i]i=l,...,JV 



"iJj=i,...,jv 



(19) 



■5 



where 7jj is as in (fT8|) , while 

1 

Vi := — 
2 



e- v ^a(x)R^ 1 (x)dx (20) 



The significance of the generating function is that the correlation functions follow by functional 
differentiation 

1 S n 

/?(„)(«!, ...,u n ) = — — -— - — — — -Z N [a] (21) 

Zpf[a\ oa{ui) ■ ■ ■ oa(u n ) a=i 

Using (fT?} and (j2"Tj) . Tracy and Widom 28J showed how the known n X n quaternion determinant 
formula from [TS] for p( n ) in the case f} = 1, N even could be reclaimed. Crucial to their method is the 
identity 

Det(l p + A pxq B qxp ) = Det(l, + B, X pA>x<,) (22) 

where A pX(? is a p x g (and B qxp is a g x p) matrix valued integral operator. However, (3=1 for N odd 
has not yielded to this approach. Yet, the correlations in this case are known — how then were they 
derived? 

2.4 Guessing the answer — then verifying it 

For general jpdfs Pn,/3, symmetric in all variables, the n-point correlation /?(„) can be expressed as an 
integral according to 



p (n) (r 1 ,...,r n ) = N(N-l)---(N-n+l) / dr n+1 ■ ■ ■ j dr N P NJ3 (r u ...,r N ) (23) 
Note in particular that 



P(n)(n, ■~,r n ) = — - - J dr n+ i p (n+ i)(ri,...,r„ + i) (24) 

This recurrence together with the initial condition 

P(N)(n, -,r N ) = Nl P Njl3 (ri, ...,r N ) (25) 

completely determines {/0( n )}n=i,2,...,iV- 

Now the functional differentiation method, when used in conjunction with a Pfaffian form of the 
generating function, is well suited for the computation of the 1- and 2- point correlations. Such a 
calculation does not rely on (|22|) . and so can be applied to both the N even and N odd cases for [3=1. 
It reveals that p/ n -\ can, for n = 1, 2 be expressed as a 2n x 2n Pfaffian with entries independent of n. 

From this, one can then conjecture a Pfaffian form for /?(„) in general. The validity of the conjecture 
requires verifying both (j2"4")l and ([23| . It is the verification of ([23)) that requires N odd be treated differently 
The verification of (|2"4")l is done by using a recursion for integrals of quaternion determinants, 
known as the Dyson Integration Theorem [7l [2ll [2] , and the parity of N plays no explicit role. 

Theorem 2.3. Dyson Integration Theorem 

Let f(x,y) be a function of real, complex or quaternion variables where 

f(x,y) = f(y,x) (26) 

with f being the function f , the complex conjugate of f or the dual of f depending on whether x and y 
are real, complex or quaternion respectively. 
Also let 



f(x,x)dn(x) = c (27) 
f(x,y)f(y,z)dp(y) = f(x,z) + Xf(x,z)-f(x,z)X (28) 
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for some suitable measure d\x, a constant scalar c and a constant quaternion A. 
Then for a matrix F nxn = [f(xi,Xj)] n xn ™e have 

J QDet[F nXn ]d^(x n ) = (c-n+l) QDet[ J F ( ;„_ 1 ) X ( n _ 1 )] (29) 

Remember that if f(xi,xj) is a scalar, then QDet reduces to a standard determinant. For a proof of 
Theorem [231 see Theorem 5.1.4 in [2"2]. 



The Pfafhan formula (|17p involves an arbitrary set of monic polynomials {Rj(x)}j=o,i.... For the im- 
plementation of Dyson's integration formula, in particular the validity of (|2"T|) and (j2"8|) , these polynomials 
must be chosen to have a certain skew-orthogonality property. 

Define a skew-symmetric inner product (-|-) by 

-I />00 /"OO 

:= - / dx e- v ^f(x) / dy e^giv) sgn(y - x) (30) 
and let {R n (x)} n =o,i,... be a corresponding family of monic skew orthogonal polynomials so that 

(i?2m|-R2ri+l) = — {R2n+l\R2m) — f n 5 nm (31) 
\R2m\R2n) = (R2m+l\R2n+l) = (32) 

With this specification of {Rj (x)}j=o,i... the correlations for the jpdf ([3]) have been shown, by verifying 
(|24| and (|25p . to be given by explicit n x n quaternion determinant formulae. First the case when ./V is 
even. 

Theorem 2.4. Let {R n {x)}n=o,i,... be as above, then for an eigenvalue distribution given by the nth 
order correlation function, in the case of N even, is given by 



P( n ){xi, ...,x n ) = QDet[f(xi,Xj)] it j=i t ... in (33) 



whe 



and 



f(x,y) = 



S(x,y) I(x,y) 
D(y,x) S(y,x) 



(34) 



N ^ e -v(y) , , 

S (. X >V) = E {<f>2k(x)R2k + l(y) - ®2k+l(x)R2k(y)) (35) 



k=o Tk 



N {j^ 1 e -V{x)-V{y) . > 

Dfoy) = ~ {R2k(x)R2k+i(y) - R2k+i(x)R2k(y)j (36) 



N/2-1 

I(x,y) = — (l>2k+i(x)<I>2k(y) - $2k{x)$2k+i(y)) +h{x,y) (37) 



k=0 



<f> k (x) = / /i(ff,x)i? fc (y)e- y ^ ( iy (38) 



= isgn(y-x) (39) 
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For N odd, the correlations exhibit the same functional form but with some modifications to the 
kernel (l34l) as follows. 



Theorem 2.5. 



with 



(AT-l)/2-l v{y) 

S odd (x,y) = T —\^k{x)R2k + i{.y)~^2k+i{x)R 2k {y) 



k=Q 
f(N-l)/2) 



■F(x)R N - 1 (40) 



D odd (x,y) = (i? 2fc (x)i? 2fe+1 (y)-i? 2fc+1 (a:)i? 2fe ( 2 /)) (41) 



(JV-l)/2-l 

/ odd (z )2 /) = £ - ($ 2fc +i(x)$ 2 fe(y) - $2A(z)$2fc+i(y)) + ») (42) 



+ ^ (6 JV _ 1 (a ! )f(y)-F(a!)d w _i(y)) (43) 



'(N-l)/ 

r„ := r n (n = 0, (N — 3)/2) (44) 

/•OO 

r ( jv-i)/2 := / e- v ^F(y)R N ^(y)dy (45) 
i?„(x) := i?»(x)- Rn -^ f°° e~ v ^F(y)R n (y)dy (46) 

r (JV-l)/2 Joo 

R N -x(x) := R N -x(x) (47) 
K{x) := / e- vto) /i(y,a ; )i?„(y)d 2 / (48) 



= i (49) 
Proofs of these theorems, using this notation, are found in [ID]. The original proofs are contained in 

us nam. 



A problem presents itself if the jpdf is not a symmetric function in all variables. Then the recursion 
([24]) does not hold and so the strategy leading to Theorems 12.41 and 12 . 51 is not valid. As mentioned, there 
is an alternative strategy involving the functional differentiation formula (|2ip . It however has not been 
successfully implemented in the case of (3 = 1, N odd. 

As to be revised in the next section, the GinOE, N odd case falls prey to both of these negative 
results, and so requires investigation of terra incognita. 



2.5 Ginibre Orthogonal Ensemble 

In the development of the theory of the Gaussian and circular ensembles, the eigenvalue jpdfs ([!]) and @ 
were contained in the pioneering papers of Dyson and Mehta 6 ( 23 . In contrast, it took over 25 years 
from the formulation of GinOE [14 to the determination of its eigenvalue jpdf. That this problem is 
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more complex than those solved previously can be seen from the fact that for any N there is a non-zero 
probability of having at least one real eigenvalue — the real line is populated despite having measure 
zero inside the support of the set of all eigenvalues. 

The jpdf therefore breaks up into sectors labelled by the number k of real eigenvalues (this number 
must have the same parity as the size of the matrix N because the complex eigenvalues come in complex 
conjugate pairs) and was first computed in [17] and then later in [8] and [24]. The explicit form reads 

k (N-k)/2 

P/v,fc(A,W) = — TTe" A * /2 TT e-K 2 +- 3 2 )/ 2 erfcdlm^)!^) A(A,W) (50) 

Cn * Li f=i 

where A = {A}^...^ are the real eigenvalues and W = {lUj, u)j}j = i i ... i (./v-fc)/2 are the complex eigenvalues, 
which come in conjugate pairs. A is the the product of absolute differences, i.e. A({xi}) := Yii<j \ x i~ x j\- 
In any one sector k the correlation function for m (where m < k) real eigenvalues and (n — m)/2 
complex eigenvalues in the upper half plane are specified by an integral analogous to (|23p . and thus 
recursions (one an integral over a real eigenvalue and the other an integral over a complex eigenvalue) 
analogous to (|24p hold. However, confining the correlations to a particular sector is not in keeping with 
the realities of the problem, in which the number of real eigenvalues is not known a priori. Consequently, 
the correct way to specify the correlations is to first introduce the summed-up generating function 

jV 

Z N [u,v] = ^2 * z k,{N-k)/i[u,v] (51) 

fe=0 

where * indicates that the sum is restricted to values of k with the same parity as N and 

k 



/oo y* oo ^ 

d\i ■ ■ ■ l d\ k Y[u(Xi) 
-OO J OO \ — \ 

„ f (N-k)/2 

x / dwx--- I dw, N _ k)/2 TT v{wi) P w ,fe(A,W) (52) 
Jr 2 , Jv. 2 , , , 



with wi = (xi,yi), wi := xi +iyi. 

From (|51[) we compute the correlation function for m real eigenvalues, and (n — m)/2 complex eigen- 
values in the upper half plane according to the functional differentiation formula 

1 6 n 
P(n)(x,w) = — — — — — rZ N [u,v] (53) 

Zn[U, V\ dU(Xl) ■ ■ ■ du{X m )OV[Wl) ■ ■ ■ dV[ W n- m ) u=v = l 

Expanding out the functional differentiation does not lead to an integral formula for p( n \ analogous to 
(|23p . because the RHS of such a formula would now consist of the sum of many different integrals. Thus, 
for the calculation of the correlations, the strategy of using (|2"4"]l and (|2"5"|) or their analogues, is not 
available for GinOE. 

This leaves then just the one strategy, namely to give a Pfaffian form of the summed-up generating 
function (p)Tj) . and to deduce from this (using (|2"2"|) or similar) a Pfaffian form for p/ n y It was some 
time after the discovery of (|50p before this program could be carried through. The first advance was the 
deduction of a Pfaffian formula for Zm\u,v\ [25] 

2~JV(iV+l)/4 

Zn ^ V ^ = ttJV f / ;/9 x Pf i a jA u } + Pj,kMi,j=h-,N (54) 

where 

a jA u ) = / dx u(x) / dy u(y) e~ (x +v ^ 2 pj_ 1 (x)p k -x{y) sgn(y - x) (55) 



f}jjt[v] = 2i I dw v(w)e y x eTic(V2\lm(w)\)(^p j ^i(w)p k ^i(w) - p k -i(w)pj-i (wj) (56) 
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for arbitrary monic polynomials pi(x) of degree 

To proceed further, the polynomials pi(x) must be chosen to be skew-orthogonal with respect to the 
skew inner product 

(pj-i,Pk-i) ■= (a 3 +/3k) 



satisfying 

G2j,2k = G2j-l,2k-l = , G2j~l,2k 

The requisite polynomials are identified in [11] 



=: Gj.k 

=1 

-G2k,2j-1 = 1"j-l$j,k 



(57) 
(58) 



p 2 j(x) = X 2] , p 2j+ 



l = x 



2j+l 



2jx 



2i-i 



r,-_ a = 2\/2^r(2j - 1) (59) 

With knowledge of these skew-orthogonal polynomials, it is a relatively routine task to then compute 
correlations using (|53|) . In the case of N even, through the use of (|22"]) . this leads to the computation of 
pi n ) , n arbitrary, for the real-real and complex-complex correlations in [11] , and then for mixed real and 
complex correlations in [3J. 

For N odd it is not known how to make use of (|22[) (recall the remark below that formula). Nonetheless, 
provided n — 1, 2 the functional differentiations required by (|53[) are tractable via a different calculation. 
Indeed, the authors have carried through the required computations. These low order calculations, 
together with the structure of the general n-point correlations for TV even as known from [3J immediately 
suggest the n-point correlations in the case of N odd. However, what is not immediate, given the 
inapplicability of the recurrence (|24p (or an analogue) and the lack of knowledge on how to implement 
functional differentiation systematically for N odd, is a way to verify these formulae. 

This suggests revisiting the problem of computing the correlations for (J5J in the case of N odd. Can 
we devise a strategy that reclaims Theorem 12.51 and will furthermore be applicable to the N odd case of 
GinOE? 

In fact this broader question is now of more importance than the calculation of the correlations for 
the N odd case of GinOE. Very recently Sommers and Wieczorek [STJ studied the correlation for GinOE 
through functional differentiations combined with the use of Gassmannian integrals. By the use of a 
certain additional 'artificial Grassmannian', they have solved the N odd case. Our idea is quite different; 
we will show that the N odd correlations can be obtained from the N even case by a limiting procedure. 
Moreover, our approach is generally applicable, when given a certain determinant or Pfaffian form for 
the correlations pertaining to the N even case. 

3 The taming of the odd 

To go from N even to N odd in (J3j we propose taking one of the eigenvalues off to infinity. This will be 
a useful strategy if the jpdf exhibits the factorisation 

P Nt i(xi,...,x N ) ~ /iv(xi) Pjv-i,i(x 2 , ... ,x N ) (60) 
Indeed ([3]) satisfies f|60p with 



fN ( x ) = 9^i x N-i e -v (X ) (61) 



It then follows from (l23l) that 



P(jn){n,-,r m ) ~ Nf N (r m ) p (m _ 1} (ri,...,r m _i) ( 62 ) 

ri | — »oo 



where the superscripts on the p(u\ indicate the total number of eigenvalues. 
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With the total number fixed on the LHS, the total is reduced by one, and thus odd, on the RHS. We 
can further use (|60|) . (|61|) and (|23l) to show 



and thus rewrite (|62)l to read 



P(m)( r l'-' r "») ~ P$)( r mV(> n -l)( r li-">*m-l) 
|r*i |— >oo 



(63) 



(64) 



Either way, we have that the (m — l)-point correlation for N odd is a limit of the m-point correlation 
for N even. Our task now is to show how (f6"2| can be used to deduce Theorem 12.51 from Theorem 12.41 
and for that we begin with humble row and column reduction. 

3.1 Humble row and column reduction 

As discussed above, the mth correlations of the (3 = 1 ensembles consist of quaternion determinants of 
m x m matrices with quaternion elements. 

Using Theorem 12 .41 we isolate the mth eigenvalue and write out the quaternion determinant explicitly, 
in Pfaffian form (according to (O). 



P(rn)(xi, -,X m ) = Pf 



I{xi , Xj ) S{Xi , Xj ) 
SyXjjX'iJ DyX<i)Xj) 

I{%m i %j ) S^Xjn , Xj ) 



S^Xjyi! Xi) D{Xi, *^m) 



o 



o 



S{x ra ^ x m ) 

for some fixed m. This matrix consists of 4 submatrices of sizes: 

• Top left: 2(m - 1) x 2(m - 1). 

• Top right: 2(m- 1) x 2. 

• Bottom left: 2 x 2(m - 1). 

• Bottom right: 2x2. 

Now since (Pf [X]) 2 = Dct[X] for any anti-symmetric matrix X, ([55)) becomes: 



S(x m: x m )Pf 



-S (xj,x^) (x^,Xj) 



-S(x rn ,Xi) 



I\3*mi %j ) S(_X m , ) 




1 " 







-1 



(65) 



l,...,m — 1 



l,j'= l,...,m— 1 



I yXi) Xj) S \Xi ^ Xj ) 



(66) 



i,j= l,...,m— 1 



n 



where 



D*(xi,xj) = D(xi,xj) 
S {xi , xj) 1=1 S\pCj l < 1 Xj^j ~ 



D^Xi , X m J S ( y X m , Xj ) SyXfji j Xj^DiXjrt , ) 



${Xm i *£m) 


^yXm ? 3>m ) 




■D\Xm ) 3£j ^)l{Xi , X m ) 


S{x m , X m ) 


S( y x 7n ^ x m ) 


$\Xi) X m )/(x m , 3Sj ) 









(67) 
(68) 
(69) 



The equality in (|66|) can be seen by using the Laplace expansion method for Pfaffians discussed in 
Section |2~T1 Note that ([66]) factors out p?^ as required by ([64]) . 

3.2 To infinity and beyond 

To compute p^^^ according to (|64p , we must cancel S(x m , x m ) and take x m — ► oo in the entries of the 
remaining Pfaffian. 

To reclaim Theorem 12.51 we must then have: 



Now using ([55]) , 



D (Xi , X j ) 
S (Xi , ) 
/ (x^ , Xj ) 

and (|37p we see that as x m — > oo 



N—>N—1 



N — >N—1 



D odd (x hXj ) 

K> 

S° dd (xi, Xj ) 
I odd {x hXj ) 



D{Xi j Xm) 

Si^Xi , x m ) 

$ {Xm i *^z) 



${Xm , Xm) 
I{xi j X m ) 



-V> ( ) p -V(x m ) 



rjv/2-i 



-i?Ar_ 2 (x;) -Rj\r_i(x m ) 



rjv/2-1 



-$^-2(2^) J2Ar-i(a; m ) 



^ e -v(*«) 

E 

fc=0 



flafcfi^i) 5 / e- v ^R 2k (y)dy 



1 



- v ^R 2k+1 (y)dy 



-iijv-i^m) t / e- v ^R N _ 2 (y)dy 

r N/2-l Z Joo 

N/2-1 x 

V - $ 2fc+1 (x l ) - / e- v ^R 2k (y)dy 

n-L 27^ 



fc=0 



-*a*(a*) 2 / e- v Wi? 2fe+1 (y)^ 



(70) 
(71) 
(72) 



(73) 
(74) 



(75) 
(76) 



(77) 



From these we find the equalities (|70|) . ([7T|) and (j72|) . and so Theorem 12.51 is reclaimed. 
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4 GinOE into the fold 



Here we show that GinOE is amenable to the same treatment; we claim that an analogy of the asymptotic 
factorisation formula (|64|) again holds true. To understand this we return to the jpdf (|50p and impose 
the ordering constraint 

Ai < A 2 < • ■ • < X k (78) 

One significance of this is that the normalisation CV.fc, which, prior to this ordering, is dependent on 
k through the factor k\ , is then Cjv,fe = Cjv,JV-fe, a function of N and N — k only [TTl [8j. The formula 
(|52| must correspondingly be modified in its domains of integration to account for these orderings. With 
this understood, we observe that for large x 



I^ Z k,(N-k)/2[U,V\ 



N — 1,N — k 



e x2 ' 2 x N 1 Z k _ h(N _ k)/2 [u,v\ 



(79) 



The reason for this is that in this limit the leading contribution comes from the functional derivative 
acting on the A& variable as ordered in ([75)) . But according to (|50p . and the remark made above relating 
to C~N,k, the jpdf contributing to each term in (|5ip keeps its functional form, but with N i— ► N — 1 and 
k i ► k — 1 and an overall factor of C "- 1 - N -" e-^V-i. 

This latter factor is the leading x — > oo form of p^(x), as can be seen from (|53[) with n = m = 1 and 
(|79"f itself. Use of (|75|) in (|53^) for general n gives as the analogue of (fM| for GinOE 



( x - w ) ~ P(i,o)( a! i)^ 1 i:i, n2 )(x,w) 



(80) 



where % is the correlation function for ni real and ri2 complex eigenvalues from the ensemble of 

N x N matrices. 

To proceed, we first need the appropriate even solution; for this we use that of [3]. 
Theorem 4.1. Let {pi}i=\ y ... be the skew- orthogonal polynomials of i59\) then 



for 
where 



< 1 ,n 2 )(x,w)=Pf 



K N (xi,Xj) K N {xi,wi) 
K N (w k ,Xj) K N (w k ,wi) 



Xi € R , u>i £ 



D(s,t) S(s,t) 
-S(t,s) I(s,t) 



K N {s,t) -- 
x = {xi, ...,x ni } , w = {wi, —,w n2 } 



(81) 
(82) 



k=0 



£>(m,??) = 2 Y] — q2k(p)q2k+i(v) ~ 92fe+i(M)92fe(??) 



2 X! _ q2k(^)T2k+l(r]) - q2k+l(n)T2k(?l) 



k=0 



2 — T"2fc(^)T2fc+i(f?) - r 2 fc+i(/i)r 2 A ; (?7) + e(jLt, 77) 



fc=0 



(83) 
(84) 
(85) 



and 
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qi(z) = e* 2 / 2 V /erfc(%/2|Im(z)|)p 4 (z) 



(86) 



n ( z ) = J ie?' 2 ^/ cr f c (V2|Im(z)|) Pl (z) z G R+ (8?) 
1 -§*<(*) z£K 



isgn(2! - z 2 ) z 1; z 2 e R 
otherwise 



$j(x) = I sgn{x - z) Pj (z)e- Z * /2 dz (89) 



oo 



The matrix in (|8ip consists of four blocks representing the four possible combinations of eigenvalues: 
real-real, real-complex, complex-real and complex-complex. For ease of manipulation we shift the mth 
real eigenvalue to the right-most column and the bottom row. Since we are shifting two rows and columns 
an even number of times we do not change the Pfaffian. We now have: 

K N (xi,Xj) K N (x u wi) K N (xi,x m ) 
Pf K N (wk,x ) K N (w k ,wi) K N (w k ,x m ) (go) 
K N ( ) K N (x wi) K N (x 

The sizes of the submatrices are: 

• top left: 2(m — 1) x 2(m — 1); top centre: 2(m — 1) x 2(n — m); top right: 2(m — 1) x 2. 

• centre left: 2(n — m) x 2(m — 1); centre: 2(n — m) x 2(n — m); centre right: 2(n — m) x 2. 

• bottom left: 2 x 2(m — 1); bottom centre: 2 x 2(n — m); bottom right: 2x2. 

It can be seen that this matrix is now equivalent to that in Equation (|65|) . By applying the same 
process as in Section |3] above, we arrive at the odd case. The details of the calculation involve index- 
shuffling with a liberal coating of tedium and so they are omitted. The results, however, are: 



Theorem 4.2. For N odd, h81}) and Il8ty) hold with the following modifcations 

(w-i) 1 

1 r 

q2k(p)q2k+i(v) - hk+i(lJ.)q2k(v) 



D{n,r,) = 2 



E 

k=0 



1 



(JV-1? 
2 

E 

k=0 



E 

k=0 



1 

rk 
1 

rk 



q2k (p)T2k+l (V) - 92fc+l {n)T2k {rf) + K(n, T]) 

T 2 k{n)T2k+i{ri) - f 2 fc + i(/x)f 2 fe(r?) +e{^,n) +8(n,r)) 



(91) 



(92) 



(93) 
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where 



Uz) = e 7 - 2 ' 2 crk(V2\Im(z)\) p t (z) 



ti(z) = 



0(ji,T)) 



Pi( z ) 



v l = 77 



* 2 / 2 yj evtc(V2\Im(z)\) pi(z) z e K 2 



z e 



o 77 e 

izz(&N-i(v)-*N-i(ii)) M,?7eR 
"W^-ifo) M 6 R , r? G Rj 



2v N 





V 




Pi (z) ~ 


1 

v N 


1 r 00 

2 7-00 




u(X)e 



/tx e RJ , v s 

/j,77<E R+ 



pjv_i(«) i = l,...JV-2 



u(A)e- A /2 ^-i(A) dA 



(94) 
(95) 

(96) 

(97) 

(98) 
(99) 



This is the exact form of the N odd correlation functions for GinOE. However, there are still two 
further points to appreciate. One is that there are inter-relations between the D,S,I within the various 
blocks, i.e. for the real-real correlations, complex-complex correlations, etc. For the real-real case the 
relationships are presented in they amount to differentiation and integration of S r ,r{x,y) to obtain 
D r<r (x,y) and I r ^ r (x,y) respectively. In the complex-complex case, the operation is complex conjugation 
of one of the variables and multiplication by i. The mixed cases are a combination of these operations. 
The inter-relationships are the same for both N even and N odd. 

The other point is that there exist summed-up forms of, say, the S in each block, and thus the D and 
I by the inter-relations, which are applicable for both N even and ./V odd. In the case of S r ^ r (x, x), this 
has been known since the work of Edclman et al. [9] , and has been shown to be true in general in [27] . 

We collect together the formulae bearing on these points in the Appendix. 



5 Concluding remarks 

The problem of computing correlation functions for GinOE is complicated by the eigenvalue jpdf breaking 
up into sectors according to the number of real eigenvalues. It is further obscured by the need to treat the 
N even and iV odd cases separately. This latter complication is shared by the general j3 = 1 eigenvalue 
jpdf ((3|). It has motivated us to revisit the problem of computing correlations in the N odd case for f3|, 
with our aim being to devise a method applicable to the GinOE for N odd. We explain, in Sections 
12.41 and 12. 5[ why the established method for computing correlations for ([3|) in the case N odd cannot be 
applied to GinOE. 

The approach we take is to compute the N odd correlations as limiting cases of the N even correlations 
according to (|64[) . in relation to ([3|), and (|80|) in relation to GinOE. Pfaffian forms of the appropriate size 
are obtained from these by row and column reductions of Pfaffian form implied by the LHSs. 

We remark that (|64| and (|80|) work equally well for relating the correlations for ./V odd to those for TV 
even, which can be used as a consistency check on our own workings. Further, there are other problems 
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to which our method applies. The most immediate is the partially symmetric GinOE |17[ I12j , which will 
be the topic of a future publication. 
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Appendix 



The explicit forms of the terms in the correlation kernel (|82[) for N odd are contained below. The 
subscripts r and c identify real-real, real-complex, complex-real and complex-complex factors. The 
convention used in this Appendix is x, y £ R and w, z £ R^. 

As mentioned after Theorem 14.21 there are inter-relationships between the various kernel elements as 
seen in [TTJ [3J [57] • Note the interesting fact that to obtain I from S you operate on the first variable, 
and to obtain D from S you operate on the second variable. 



v i 

I r<r {x,y) = I S r , r (z, y)dz + -sgn(x - y) 



d 

D r , r (x,y) = -—S rtr (x,y) 
ay 



Ir t c(x,w) = -I Cir (w,x)=— S TtC {z,w)dz 

J x 

D rtC (x,w) = -D c<r (w,x) = —iS rc (x,w) 



(A.l) 

(A.2) 

(A.3) 
(A.4) 



I Ct r(w,x) = -I riC (x,w) = iS c ,r{w,X) 

d 

D c , r (w,x) = -D riC (x,W) = - — S c ,r(w,X) 



I CyC (w,x) = iS c ^ c (w,z) 



dx 



(A.5) 
(A.6) 
(A.7) 



D CiC (w,z) 



-iS c>c (w, z) 



(A. 



Below are listed the kernel elements S for real-real, real-complex, complex-real and complex-complex 
correlations. Using the above relationships, the other elements can be deduced. 



S r ,r(x, y) 



72 



E- P2k+l(x)$2k(y) ~ P2k(x)^2k+l(y) 



k=0 



Tk 



(p2k+l(x)$N-l(y) - PN-l(x)$2k+l(y)j 

+ 2fc+2 (p2k(x)$N-l(y) - PN-l(x)$2k{y) ) 
VN V / 



V2k+l 
VN 

V N 
-x 2 /2 



2v N 



-p N -i(x) 



S r<c (x,w) = 2ie- x2/2 e- a2 / 2 ^erfc(V2\lm(w)\) 



N-l 



-1 



X — P2k{x)p2k+l{w) ~ P2k+l(x)p2k(w) 



k=0 
V2k+1 
1> N 
'2k^ 



(pN ~l{x)p2k+l{w) ~ P2k+l{x)pN-l{w) 

+ V2k+2 / _ 1 (x)p 2 k(w) ~ P2k{x)PN-l{w) J 
VN V / 



(A.9) 



(A.10) 
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S c . r (w, X) 



-w 2 /2 



erfc(v / 2|Im( u; )l) 



N-l 



-1 



X ^ — P2k+l(w)$2k(x) - P2k(w)$2k+l(x) 



k=0 
V2k+1 

V N 
V2k+2 

V N 
p -w 2 /2 



Jerfc(v^|Im(u;)|) p N -i(w) 
S c , c (w,z) = 2ie-" ,2/2 e-* 2 / 2 -\/erfc(v / 2|Im(w;)|)Y / erfc(v / 2|Im(z)|) 

J ^ r- 

X V" — P2fe(w)p2fc+l(z) -P2/c+l(^)P2fc(2) 
fe=0 



^2fe+l 
1> N 



p N ^l(w)p 2 k+l(z) - P2k+l{w)p N -l(z) 



+ _2fc+2 / N _ x {w)p2k{z) ~ P2k{w)pN~l{z)) 
V N \ / 



(A.11) 



(A.12) 



The kernel elements S can also be summed-up |12l I27j ; the summed up forms are listed here. Note 
that these formulae are insensitive to the parity of N (in fact N need not even be an integer), and so 
provide a further check of our work in this paper. 



s r , r ( x >y) 

S r ,c{x,w) 



1 r -(x-yf/2 T(N-l,xy) _ (JV -3)/2 -x 2 /2 N-l 7( V: V /2) 

r(iv-i) ' r(A^-i) 



' 27T 

ie -(x-w) 2 /2 



1 r -( w - x) 2 /2 r(iv- i,wx) _ 2 (at_ 3 )/2 p -^ 2 /2 /v-i 7(^-^ 2 /2) 
r(JV-l) ' ' T(N-l) 



erfc(v / 2|Im(w)|) 



(A.13) 
(A.14) 



(A.15) 



S c ,c(w, z) 



' /2 < 2 ~ w ^ { t (N ZYT \feMV2\Im(w)\)^/eMV2\Im(z)\) (A.16) 



2tt 
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